Partial order and finite-temperature phase transitions in Potts models on irregular 

lattices 
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FIG. 1: (Color online) (a) Union- Jack lattice. Sites in sub- 
lattices A (red circles) and B (blue) have coordination num- 
bers za — zb — 8, while those in sublattice C (yellow) have 
zc = 4. (b) Diced lattice. A sites (red) have za = 6, while B 
sites (blue) have zb = 3. (c) Centered diced lattice: za = 12 
(red), zb = 6 (blue), and zc = 4 (yellow). 
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We evaluate the thermodynamic properties of the 4-state antiferromagnetic Potts model on the 
Union-Jack lattice using tensor-based numerical methods. We present strong evidence for a previ- 
ously unknown, "entropy-driven," finite-temperature phase transition to a partially ordered state. 
From the thermodynamics of Potts models on the diced and centered diced lattices, we propose that 
finite-temperature transitions and partially ordered states are ubiquitous on irregular lattices. 
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The g-state Potts model has an essential role in the 
theory of classical critical phenomena and phase transi- 
tions jlj. Antiferromagnetic Potts models are rich and 
complex, displaying many different types of behavior as 
a function of q and of lattice geometry, and the search 
for guiding principles continues. Key questions include 
whether the model has a phase transition, if this occurs 
at finite temperature, the nature of the low-temperature 
phase, and the universality class of the transition. 

In the Landau theory of phase transitions, these are 
described by an order parameter. Minimizing the en- 
ergy at low temperatures requires the order parameter 
to be finite everywhere. However, in systems with ex- 
tensive ground-state degeneracy [||, of which the best- 
known example is ice the nonzero entropy at zero 
temperature may cause a different type of ordered state. 
A partial order, involving only some of the lattice sites, is 
stabilized by minimization of energy in combination with 
maximization of entropy. This "entropy-driven" transi- 
tion can occur at a finite temperature. Such partial order 
is known in both frustrated [Q, || and unfrustrated || sys- 
tems. In the latter case, partial order arises purely due 
to configurational entropy effects, and the q = 3 antifer- 
romagnetic Potts model on the diced lattice provides an 
excellent example of the associated phase transition 0. 

In this Letter we pursue the physical origin of the 
finite-temperature phase transition in two-dimensional 
(2D) q-state Potts models. Although these systems have 
no exact solution for q > 2, we employ tensor-based nu- 
merical methods to obtain hitherto unavailable thermo- 
dynamic quantities. We show that the q = 4 Potts model 
on the Union- Jack lattice exhibits a finite-temperature 
transition to a state of partial order. We characterize 
this transition by computing the entropy, specific heat, 
and magnetization, quantities we also use to provide a 
complete discussion of the q — 3 model on the diced lat- 
tice. We propose that finite-temperature transitions and 
partially ordered states are a general property of Potts 
models on irregular lattices. 

Let the q states of the Potts model for lattice site i be 



labeled Ui = 0, 1, q — 1. In the Hamiltonian, 



(1) 



J > corresponds to the antiferromagnetic case and a 
field H is coupled to one of the q states. We consider only 
a single coupling J on every bond, and begin with H = 0. 
We employ tensor-based numerical techniques developed 
recently by a number of authors [p|-ft2| to compute the 
partition function Z to high accuracy, and hence to ob- 
tain the required thermodynamic quantities. 

We illustrate the methods by focusing on the Union- 
Jack lattice [Fig. 1(a)]. The partition function 



Z = Tre-^ = J2Ii e ~ mD > 
□ represents the structural unit (Fig. 2(a)) and 



(2) 



(3) 



The tensor T n = e~ m ^ is denned ||, [T| by sum- 
ming over the C-sublattice sites (T5 and introducing the 
bond (or dual-lattice) variables a = mod(cri — 02, 9), 
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FIG. 2: (a) Schematic representation of the tensor T, ob- 
tained by integrating over the C-sublattice sites and per- 
forming the dual transformation, (b) The Union- Jack lattice 
(dashed lines) is transformed to the dual square lattice (solid 
lines) on which the tensor T is defined. 

j3 = mod(<72 — (73,(7), 7 = mod(<73 — a 4, q), and S = 
mod((74 — ai,q) [Fig. 2(a)]. The partition function 



a/37<5. . . 



(4) 



becomes a product of tensors defined on the square lat- 
tice [Fig. 2(b)]. The tensor renormalization group (TRG) 
H [U)[ is a real-space coarse-graining method, in which 
the precision of the tensor contraction is greatly enhanced 
by simultaneous renormalization of an "environment" 
block 10, O . Infinite time-evolving block-decimation 
(iTEBD) Jl2| is a projection method, in which projecting 
the transfer matrix sufficiently many times on a random 
vector gives very accurately its largest eigenvalue. Both 
methods simulate the properties of an infinite system, the 
truncation being performed in the size D of the rank-4 
tensor T a j3 7 s, which thus determines the precision. With 
an accurate partition function, we then obtain all other 
thermodynamic information. 

One guiding principle for the antiferromagnetic Potts 
model in 2D concerns the existence of a critical q, q c (C), 
for each lattice geometry C. At zero temperature, neigh- 
boring sites may not have the same "color" oi (Fig. 1), 
making the model equivalent to a vertex coloring prob- 
lem. By exploiting the Dobrushin Uniqueness Theorem 
p3[ , Salas and Sokal proved [Q that the correlation func- 
tion decays exponentially at all temperatures (includ- 
ing zero) for sufficiently large q, meaning a disordered 
ground state and no phase transition. For sufficiently 
small q, ordered ground states usually exist |?], fl5j. For 
q = q c (C), T = is a critical point; this situation is com- 
mon for many integer q values, including on the square 
p5| , kagomc jl6| (both q c = 3), and triangular (q c = 4) 
lattices £l7|. 

In terms of q c , the diced lattice [Fig. 1(b)] is anoma- 
lous. Despite an average coordination number of 4, the 
model with q = 3 has a finite-temperature phase transi- 
tion, hence q c > 3. This transition is driven J?J by the 
entropy available from the irregular nature of the lattice, 
meaning that there are more sites of one type than of oth- 
ers. The low-temperature phase can be ordered on the 
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FIG. 3: (color online) (a) Entropy and specific heat for the 
q — 4 Potts model on the Union-Jack lattice, computed with 
D — 800. Thin, green lines denote the zero- and infinite- 
temperature limits of the entropy. Inset: specific-heat dis- 
continuity near T c . (b) Linear fit of T c with 1/D, which ex- 
trapolates to T c (oo) = 0.339(1). The error bar denotes the 
upper and lower temperatures at the discontinuity. 



A sublattice (for example, Ui = 0) but not on B (<x; = 1 
or 2), creating a partially ordered state. We investigate 
the hypothesis that such behavior is generic for irregular 
lattices. We do this both by seeking an explicit new ex- 
ample of a finitc-tcmpcraturc transition, which leads us 
to consider the Union- Jack lattice, and by analyzing in 
detail the thermodynamic properties of irregular lattices. 

We begin by considering the q = 4 Potts model on the 
Union- Jack lattice. This centered square lattice is tri- 
partite [Fig. 1(a)]. While the q = 2 (Ising) model shows 
partial order and the q — 3 model is perfectly ordered, 
the q = 4 model may be expected from its average co- 
ordination (z = 6) to have zero-temperature order. To 
address the question of a phase transition occurring in- 
stead at finite temperature due the irregular nature of 
the lattice, we focus directly on the specific heat, shown 
in Fig. 3(a). The small but clear gap in the curve [in- 
set, Fig. 3(a)] indicates a discontinuous second deriva- 
tive of the free energy, and hence a second-order phase 
transition. We have performed detailed calculations to 
delineate the nature of the discontinuity over a range of 
D values, and exploit the linear scaling behavior of the 
matrix-product state with 1/D |ll|, shown in Fig. 3(b), 
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to obtain the precise result T c — 0.339(1). Thus the 
g = 4 model on the Union- Jack lattice does indeed pos- 
sess a finitc-tcmperature phase transition, becoming only 
the second known Potts model in this category. 

This transition was to our knowledge neither known 
beforehand nor even predicted. It underlines directly 
the importance of the irregular nature of a 2D lattice 
in promoting a finite zero-temperature entropy and par- 
tial order in the ground state. Our result is connected 
to several other models in statistical mechanics. First, 
some Potts models may be mapped to a height model, 
and when this mapping exists the system is critical or 
ordered at zero temperature (?], [ll|. As expected from 
our result, such a height map does exist for the q = 4 
Union- Jack lattice jl7|]. Second, at zero temperature, 
the q — 4 Potts model on the Union- Jack lattice can be 
mapped directly to the 3-bond coloring problem on its 
dual, the 4-8 lattice |l7|]. The total number of states 
on the 4-8 lattice is known and hence the zero- 

temperature entropy of the Union- Jack lattice should be 
5 UJ (0) = 21nW 4 - 8 = 0.430997, with W 4 - 8 = 1.24048. 
Our numerical result is 0.430999 [Fig. 3(a)]. 

Third, the bond-coloring problem on the 4-8 lattice is 
further equivalent to a fully-packed loop (FPL) model 
(on the same lattice), obtained by considering all con- 
figurations of noncrossing closed loops. The FPL parti- 
tion function is Z = ^2 G tl Nl , where n is a loop fugacity 
(weight), Nl is the number of loops, and G denotes all 
loop configurations. Like the FPL model on the square 
and honeycomb lattices, n = 2 for the 4-8 lattice, but un- 
like these cases the 4-8 lattice is not critical p0| . This is 
again consistent with our demonstration that long-range 
order is present at finite temperatures. Finally, the 3- 
bond coloring model on the 4-8 lattice is also equivalent 
to a 3-vertex coloring model on the square-kagome lat- 
tice [0, and thus the q = 3 Potts model on the square- 
kagome lattice is equivalent at zero temperature to the 
q = 4 Potts model on the Union- Jack lattice. 

For a full understanding of the finite-temperature tran- 
sition and the partially ordered ground state, we turn to 
the thermodynamic quantities extracted from the parti- 
tion function. For a heuristic understanding of the par- 
tially ordered state, we begin with the entropy. The q = 3 
Potts model on the diced lattice [Fig. 1(b)] is the proto- 
typical model in this class for a ground state with partial 
order (Rcf. Q and references therein). As above, if the 
minority sites order, the majority sites would have two 
remaining degrees of freedom (d.o.f.s), giving an entropy 
per site 5^(0) = 21n2/3 = 0.462098. The entropy we 
compute is shown in Fig. 4, and its low-temperature limit 
is S D (0) = 0.473839. This minor deviation indicates that 
an ideal A-sublattice order is rather close to the true 
ground state, with only small contributions from states 
of imperfect A order. This partial order is destroyed at 
T c by thermal fluctuations, and at high T, where all sites 
may explore all three d.o.f.s, the entropy approaches In 3. 




FIG. 4: (color online) Solid lines: entropy S(T) and specific 
heat C(T) for the q = 3 Potts model on the diced lattice, 
calculated with D — 100. The dotted line marks the result of 
Ref. j?J]. Thin, green lines mark the low- and high-T entropy 
bounds. Dashed lines: results for the q = 4 Potts model on 
the centered diced lattice. 

Analogous considerations for the Union- Jack lattice 
give again a maximum entropy if the high-coordination 
(minority) sites order. This order may involve sublat- 
tices A, B, or both A and B simultaneously. In the last 
case, if a a = and ob = 1, the d.o.f. in uc = 2,3 
yields a total of 2 N / 2 states. If only A sites are or- 
dered {a a — 0, ctb,<tc = 1, 2, 3), the B and C sublattices 
form a q = 3 Potts model on the decorated square lat- 
tice. Let the zero-temperature entropy of this model be 
S D (0) = ln£, then the number of states with partial or- 
der on A only is £ 3Ar / 4 , and these will dominate the total 
if C > 2 2 / 3 = 1.587401. We have calculated separately 
the value 5 DSL (0) = 0.561070, whence £ = 1.752547. 
Thus the ground state is indeed composed primarily of 
states with one ordered sublattice. The discrepancy be- 
tween 3S DSL (0)/4 = 0.420802 and our result [Fig. 3(a)], 
S' UJ (0) = 0.430999, can be ascribed to states with neither 
A nor B order. The high-T entropy approaches In 4. 

We now return to the specific heat and to the phase 
transition. Results for the q = 3 model on the diced 
lattice are shown in Fig. 4. The transition point may be 
obtained with high precision and we find T c = 0.505(1), 
a value consistent with the alternative approach used in 
Rcf. [Q . The peak feature is much more pronounced than 
in Fig. 3. This reflects the relative lack of competition 
between different types of partially ordered state in the 
diced lattice (A order only) as compared to the Union- 
Jack lattice (A or B order). We note that both transitions 
occur near the peak of the specific-heat curve at T ~ J/2, 
which is the characteristic energy scale of the system, 
and we suggest that this behavior is generic: a finite- 
temperature transition must occur at a value T c /J ~ 
O(l), and cannot occur arbitrarily close to T = 0. 

The other thermodynamic quantities we illustrate here 
are the magnetization and susceptibility. By consider- 
ing finite fields H in the Hamiltonian of Eq. ([!]), we de- 
duce the spontaneous magnetization on sublattice A of 
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FIG. 5: (color online) Magnetization M(T) of sublattice A for 
the q = 4 Potts model on the Union- Jack lattice, calculated 
for an applied field of 0.000025. Thin green lines indicate the 
low- and high-T bounds. Inset: magnetic susceptibility x(T). 



the Union- Jack lattice [red circles in Fig. 1(a)], M 
J2ieA &<ri=o/N, from the expression 



M = -^- lim — InZ, 

dH N^oo N 



(5) 



computed with a very small external field (H = 
0.000025). The results (Fig. 5) show a clear step around 
T c . The high-T limit of M for a completely disordered 



state is M UJ (oo) = \ X \ = ^. The low-T limit for a 



state of perfect A-sublattice order would be Mq 3 (0) = j , 
and in fact this value must be obtained for temperatures 
below the energy from the applied field. However, it is 
clear from our entropy calculation that the ground state 
is not one of perfect A order. We find that the mag- 
netization exhibits a related zero-temperature deviation, 
tending towards the value M UJ (0) = 0.2232. This de- 
viation is somewhat smaller for the diced lattice, where 
the ideal magnetization would be ^ and our numerical 
result is 0.3192. Returning to the transition, our cal- 
culations show the same behavior as in Fig. 3(b), that 
T c falls slowly with increasing D. The best qualitative 
indication of the transition is provided by the magnetic 
susceptibility, \ = dM/dH, which has a robust peak 
(inset, Fig. 5). The critical exponents of these quantities 
may also be computed by the same techniques, but the 
highly numerically demanding task of obtaining adequate 
precision remains in progress. 

When the discussion of 2D antiferromagnetic Potts 
models is framed in terms of q c {£), the diced lattice is 
the only known system with an average z of 4 but q c > 3. 
The Union- Jack lattice is now revealed as the only sys- 
tem known with an average z of 6 but q c > 4. In fact this 
is the largest value known for any planar lattice. Phase 
transitions in different Potts models belong in general 
to different universality classes. While the universality 
class for the q = 4 Potts model on the Union-Jack lat- 
tice remains under investigation, the partial order we find 



breaks both the 4-state (Potts) symmetry and the (Ising- 
likc) sublattice symmetry between A and B. 

We end this discussion with the logical extension of 
our analysis. The 11 regular lattices (all sites equivalent) 
obtained by tiling the plane with regular polygons are 
known as Archimedean. Among their dual lattices, three 
are regular and the other eight are irregular. This is the 
set of Laves lattices, which includes the diced and Union- 
Jack lattices. We propose that for each Laves lattice with 
an integral average coordination number, there exists a 
Potts model with integral q that would feature a zero- 
temperature transition on a regular lattice of the same 
coordination, but has a finite-temperature transition, to 
a state of partial order, on this irregular lattice. As an 
example we cite the -D(4, 6, 12) lattice, or centered diced 
lattice, shown in Fig. 1(c). This tripartite lattice has an 
average coordination 2 = 6. Indeed we find (Fig. 4) that a 
g = 4 Potts model on this lattice has a very robust finite- 
temperature phase transition to a state of predominantly 
A-sublattice order. Demonstrating the existence of the 
same physics on a third lattice in this class very strongly 
reinforces our proposal. 

To conclude, we have demonstrated the existence of a 
previously unknown, finite-temperature phase transition 
in the q = 4 Potts model on the Union- Jack lattice. This 
establishes the essential property that the presence of 
inequivalent sites, leading to a nontrivial entropy, drives 
the finite-temperature transition and confers unusually 
high values on q c . We find this type of transition on 
other irregular lattices in 2D. Our analysis underlines the 
utility of tensor-based numerical methods in investigating 
the physics of classical statistical mechanical models. 
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